Electron channels in biomolecular nanowires 
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We report a first-principle study of the electronic and conduction properties of a quadruple- 
helix guanine wire (G4-wire), a DNA-derivative, with inner potassium ions. The analysis of the 
electronic structure highlights the presence of energy manifolds that are equivalent to the bands 
of (semi) conducting materials, and reveals the formation of extended electron channels available 
for charge transport along the wire. The specific metal-nucleobase interactions affect the electronic 
properties at the Fermi level, leading the wire to behave as an intrinsically p-doped system. 
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I. INTRODUCTION 

Fueled by the ever increasing drive for miniaturiza- 
tion and improved performance in electronic devices 
and by potential applications in the nanotechnologies, 
the research effort to investigate the properties of novel 
nanowire materials is undergoing an impressive growth. 
It is foreseen that the conventional solid-state technol- 
ogy could be replaced by new generations of devices 
based on molecular components, which take advantage 
of the quantum mechanical effects that rule the nanome- 
ter scale. Molecular electronics^ is currently explored as 
a long-term alternative for increasing the device density 
in integrated circuits. However, to keep up with these 
expectations, the new trend must provide flexible, repro- 
ducible and well structured architectures, easy to wire in 
a programmable way. 

By virtue of their recognition and self-assembling prop- 
erties, DNA molecules seem particularly suitable to ful- 
fill these requirements. Both the intrinsic combinatory 
principles of nucleic acids and their chemistry can be ex- 
ploited to build precise, miniaturized and locally mod- 
ulated patterns, where the drawing of functional arrays 
is obtained through a series of programmed chemical re- 
actions and not by the physical handling of the sam- 
ples. The realization of DNA-based wired architectures 
via self-assembly is a viable route to scale down the size 
of devices to the molecular level However, whereas 
it was demonstrated that the self-assembling capabil- 
ities of DNA make it suitable as a template to wire 
metallic materials^ its ability as an intrinsic conduc- 
tor is questioned by experimental results^ Depending 
on base sequence, molecule length, environmental condi- 
tions, substrates, and electrode materials, the direct mea- 
surements of the dc conductivity of DNA-based struc- 
tures in solid-state devices^ report insulating character^ 
semiconductor-like transport characteristics^ ohmic be- 
havior, 10 and proximity-induced superconductivity^ In- 
deed, even in the cases in which charge transport has 
been observed, the current is very low, with resistances 
of the order of 0.1-1 Gil cross the length of the DNA 



molecules (variable between 10 /zm and 10 nm). Thus, 
due to its apparent poor intrinsic conductivity, DNA 
might be reasonably considered as a bad insulator rather 
than a viable electrical molecular wire. 

Besides the standard DNA, other nucleotide-based he- 
lical molecules, such as guanine quadruple helices (G4- 
wires) or "metal- manipulated" duplexes, may offer the 
desired mechanical, recognition, and self-assembly prop- 
erties, that make DNA so attractive. With respect to 
native DNA, these derivatives have metal cations in the 
inner core of the base stack. Whereas the interactions be- 
tween external ions and the double helix have been largely 
studied both experimentally and theoretically, the effects 
of their inclusion inside the helix are largely unknown. 
The presence of internal metal ions may drastically af- 
fect the bonding pattern with and among the bases, in- 
troducing novel features in the structural and electronic 
properties of the system^ A promising pathway for the 
exploitation of DNA as a conductor in molecular devices 
is indicated by the evidence that metal ions incorporated 
in the helix core may modify the conductivity of DNA- 
based wiresiiiii4 

In this paper, we focus on nanowires known as G4- 
wiresiii (or quadruplexes) , consisting of stacked guanine 
(G) tetrads (G4): the structure of these systems is illus- 
trated and described in Figure 1. These one-dimensional 
polymers are becoming appealing as prospective candi- 
dates for bio-molecular electronics because, due to the 
low ionization potential of guanine (the lowest among 
nucleic-acid bases), they might be suitable to medi- 
ate charge transport by hole conduction along the he- 
lix, and have even been suggested as nano-mechanical 
extension-contraction machines^ In the presence of ap- 
propriate metal cations (especially K + and Na + ), solu- 
tions of homoguanylic strands in water jiLi£ as well as 
lipophilic guanosine monomers in organic solvents^ self- 
assemble in right-handed quadruple helices. The G4- 
motif has been identified in both cases; however, while 
the quadruplexes obtained from guanylate strands have 
an outer mantle of sugars and phosphates (as in DNA) 
that connects the adjacent G4 planes, those obtained 
from lipophilic guanosine derivatives have no interpla- 
nar connection: a continuous linker between consecutive 
tetramers is not necessary to form the wires. The guanine 
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FIG. 1: (a) Side view of a short quadruplex from which 
the simulated G4-wire was extracted. The image was derived 
from an experimental X-ray structure.—- The central core of 
the helix (red) is constituted of stacked supramolecular struc- 
tures, the G-quartets (G4), and the K + ions (blue spheres) are 
hosted in the central cavity of the G-quartet stack. The lateral 
gray sticks describe the external sugar-phosphate backbone, 
excluded in our simulations. By virtue of the overall C4 sym- 
metry of the helix, it is possible to extract a periodic unit of 
three G4 tetrads from this finite quadruplex, that constitutes 
the elementary building-block of the simulated infinite wires, 
(b) Isolated G-quartet that forms one plane of the stacked G4- 
wire. Each G4 consists of four co-planar H-bonded guanines, 
arranged in a square-like configuration. The image shows the 
relaxed geometry and an isosurface plot of the HOMO mani- 
fold, (c) Scheme of the bypiramidal K + -0 coordination. 



quadruplexes have recently attracted interest because of 
their possible role in biological systems*^ their biological 
relevance has propelled a large number of investigations 
(e.g. X-raji2i and NMR-2) aimed at characterizing such 
unusual supramolecular structures. Quadruple helices 
have been obtained in the presence of different monova- 
lent (K+, Na+, NH4+)±2 a ^ and divalent (Ca 2+ , Ba 2+ , 
Sr 2+ , pb 2 +)iS h iSi cations. Despite the different chemical 
nature of the constituents, all these G4- wires are charac- 
terized by an inner core of stacked G4's intercalated by 
metal cations: one cation in every tetrad and one cation 
in every other tetrad, in the case of mono- and di-valent 
ions, respectively. The self-assembly capability of the G4- 
helices allows for the formation of quite long (100-1000 
nm) and stiff wires. AFM images 2 ^ of four-stranded he- 
lices obtained from G-rich oligonucleotides on substrates 
revealed the G4-wires to be uniform and relatively rigid 
polymers, with few bends, kinks, or branches. 

Despite the large amount of structural investigations, 
the conduction properties of these nanowires are basically 
unknown and a direct measurement of electrical proper- 
ties of G4-wires is still missing. From the theoretical 
point of view, quantum chemistry and molecular dynam- 
ics studies focused on the energetics and on the geometry 
of isolated G-quartets or finite clusters of stacked G4's,- 2 ? 



whereas the electronic properties of these materials were 
so far investigated to a much lesser extent £L 

In the following, we present a first-principle investiga- 
tion of the electronic and conduction properties of pe- 
riodically repeated G4- wires. Recently*^ we described 
the structure and the energetics, as well as some basic 
features of the electronic structure, of an infinite G4- 
wire with and without the presence of K + ions in the 
inner cavity. In this paper, we focus entirely on the elec- 
tronic properties of the same system and present a com- 
plete and thorough analysis of the guanine- guanine and 
metal-guaninc interactions, discussing the effects on the 
conduction properties of the tubular system. 



II. COMPUTATIONAL APPROACH 

We performed ab-initio calculations of the electronic 
and conduction properties of infinite G4- wires in the pres- 
ence of K + ions, within the Density Functional Theory 
(DFT) approach, 28 using the PW91— gradient-corrected 
exchange-correlation functional. DFT is lately gaining a 
large credit in the scientific community as a reliable and 
accurate method to describe large-scale biomolecular ag- 
gregates^ including guanine-based stacks^ 

Our total-energy-and-force calculations^ allowed us to 
attain a simultaneous description of the optimized atomic 
configuration and of the corresponding electronic struc- 
ture for the selected systems (see Figure 1): the finite 
planar guanine tetrad (G-quartet) and the infinite heli- 
cal G4-wire filled with K + ions in the inner cavity (label 
3G4/K+). The single-particle electron wave- functions 
were expanded in a plane-wave basis set with a kinetic- 
energy cutoff of 25 Ry. Two special k-points in the ir- 
reducible wedge were employed for Brillouin Zone (BZ) 
sums in the case of the G4-wire. The infinite helix was 
simulated by repeated supercell containing three stacked 
G4 tetrads^ employing periodic boundary conditions in 
the three spatial direction a thick vacuum layer (~ 16 
A) in the directions perpendicular to the helical axis pre- 
vented spurious interactions between adjacent replicas of 
the wire. For the isolated G-quartet, the same vacuum 
thickness was employed also in the third direction per- 
perdicular to the plane of the tetrad, and only the T point 
was used in the BZ sampling. 

The electron-ion interaction was described by non- 
norm-conserving pseudopotentials 33 for all the species 
(C, N, O, H) except K, for which a norm-conserving 
pseudopotential 3 ^ was used. For the latter species, both 
the valence 4s and the semi-core Zp shells contributed 
to the system with valence electrons. This treatment 
represented a significant refinement towards a complete 
description of the electronic structure, with respect to 
the simplified results^ obtained by fixing the 3p shell in 
the frozen core and by applying Non Linear Core Cor- 
rection^- (NLCC) to account for partial core relaxation. 

The starting atomic configuration was obtained from 
the results of the X-ray analysis^ 1 of the d(TG4T) 



quadruple helix. Motivated by the observation that G4- 
wires form withiLi 8 . and without^ the covalent skeleton 
and by our specific interest in the base stack as a channel 
for charge mobility, we neglected the external backbone 
in our simulations and focused on the central core of the 
helix, constituted of guanines and metal cations (see Fig. 
la). This choice is supported by theoretical reports^ 
for approximated DNA structural models, which assert 
that if any current flows in such systems, it does through 
the base-stacking, without involving the external man- 
tle in the transport phenomena. The same evidence was 
more recently confirmed by DFT calculations of the elec- 
tronic properties of real DNA sequences (A-DNA 8 . and 
Z-DNA 37 ): these simulations showed that the orbitals 
related to the sugar-phosphate backbone are a few eV 
below the Highest Occupied Molecular Orbital (HOMO) 
and above the Lowest Unoccupied Molecular Orbital 
(LUMO) of the system. In principle, one should expect 
that the environment surrounding the base stack plays an 
important role in the overall conductivity properties of 
DNA molecules! 37 ! 38 ! 39 Recent studies^ 9 . 6 pointed out 
that a random distribution of counter-ions in the unit 
cell modifies the electrical properties of DNA both re- 
ducing the bulk quantum conductance and introducing 
localized empty states in the energy gap. However, this 
result is coherent with a picture where the external ions 
quantitatively affect the global transport properties of the 
system, but do not constitute an alternative pathway for 
the electron/hole transport through the helixii There- 
fore, we believe that the bare guanine core will be well 
representative of the essential electronic properties of the 
quadruplexes, allowing for the inspection of the charge- 
migration mechanisms and of the key features of metal- 
molecule interaction. An explicit account for the effects 
of the backbone and of the surrounding counterions 3 ^^ 9 . 
would instead be demanded for a quantitative evaluation 
of the quantum conductance to be compared with mea- 
sured transport characteristics, which is way beyond the 
purpose of the present work. 



III. RESULTS AND DISCUSSIONS 

By means of the first-principle approach outlined 
above, we optimized both the isolated G-quartet and the 
periodic G4-stack. 

The analysis of the electronic structure of the planar 
G-quartet shows interesting features that will help un- 
derstanding the electrical properties of G4- wires. As also 
found for the other planar G-aggregates (e.g. dimers, rib- 
bons)^ the H-bonds among the guanines do not favor 
the formation of supramolecular orbitals extended on the 
whole G-quartet and the existence of dispersive bands. 
On the other hand, the intermolccular interactions split 
each guanine energy level into a multiplet structure: each 
multiplet is composed of four (the number of G's in the 
tetrad) energy levels and has a total width of about 200 
meV. The orbitals that contribute to a manifold have 
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FIG. 2: (a) Band-structure along the GA direction parallel 
to the axis of the 3G4/K + G4-wire, which coincides with the 
direction along which the planar tetrads stack, (b) Total and 
Projected Density of States (DOS) of the wire. The shaded 
gray area represents the total DOS. The thin red (thick blue) 
solid curve marks the projection of the DOS onto G (K). Note 
that the total DOS and the G-DOS coincide in a large part of 
the spectrum, where the thin red curve almost coincides with 
the border of the shaded area. The dashed horizontal curve 
indicates the Fermi level. 



identical character and are localized on the individual 
G molecules. Figure lb shows an isosurface plot of the 
convolution of four electron states (the HOMO's of the 
four guanines in the tetrad) which form the 7r-like HOMO 
manifold. 

By exploiting the square-symmetry of the planar 
tetrad and the 30° twist angle between consecutive 
tetrads, we simulated a quadruplex of infinite length with 
a periodically repeated unit supercell containing three 
stacked G-quartets and three intercalated K + ions in 
the unit cell^£ Each potassium ion was symmetrically 
located between two consecutive tetrads, and bipyrami- 
dally coordinated with the eight (four above and four 
below) nearest-neighboring oxygen atoms (Fig. lc). The 
atomic positions were relaxed until the forces vanished, 
within an accuracy of 0.03 eV/A. The structure, the en- 
ergetics, and the metal-induced stability of the tube were 
described elsewhere^ We find now that the explicit in- 
clusion of the semi-core 3p electrons of K in the valence 
shell for the pseudopotential calculations does not alter 
the results obtained previously within the NLCC approx- 
imation^ and does not change the understanding of the 
system from the structural point of view. The refined 
treatment of the semi-core 3p electrons of K allows us to 
gain a deeper insight into the electronic structure. 

The band-structure of the G4-wire 3G4/K + is shown 
in Figure 2a: the unoccupied states (above 3.5 eV) are 
separated from the occupied states (below eV) by a 
large energy gap (affected by the typical DFT underes- 



tiruationas), and some of the states around zero energy 
are partially filled. We previously demonstrated^ that 
7T-7T coupling among guanines may give rise to delocal- 
ized Bloch-type orbitals, whose band dispersion along 
the stacking direction depends on the relative rotation 
angle between nucleobases in adjacent planes. Taking 
into account the detected rotation angle^i of 30° for the 
guanine quadruple helices, we now report that the elec- 
tronic band-structure of the 3G4/K + column is repre- 
sented by dispersionless energy bands, indicating that no 
complete supramolecular orbital derealization is to be 
expected for such a degree of helicity. This behavior leads 
to the conclusion that in G4-wires the ir-ir superposition 
is not sufficient to induce a sizable coupling between the 
molecular states of guanines and/or the atomic orbitals 
of potassium, which would lead to purely dispersive bands 
similar to those of (semi) conducting materials. However, 
as for the isolated G-quartet, the band-structure of the 
3G4/K + column identifies the presence of manifolds that 
are equivalent to effective bands: each manifold contains 
a number (or multiplct) of levels that can be explained 
in terms of the number of nucleobase molecules and of 
potassium ions in the periodic supercell, as we discuss 
later. The orbitals in a multiplet have identical symme- 
tries and some degree of hybridization between neigh- 
boring bases. The levels within a manifold are so close 
(dense) in energy (AE < 20meV) that an external weak 
interaction (e.g. the thermal fluctuations or an exterinal 
electric field) may let the orbitals mix, giving an effective 
dispersive band. This interpretation suggests a descrip- 
tion of the electronic structure where the spreading of the 
energy levels of the manifold leads to the formation of 
effective band-like peaks in the Density of States (DOS) 
(Fig. 2b, shaded gray area). The DOS of the 3G4/K+ 
wire thus appears similar to those of materials with a 
dispersive band-structure. 

It is important to note that there is a one-to-one cor- 
respondence between the band structure and the ground 
state conduction properties (e.g., the quantum conduc- 
tance spectrum) of a periodic system in the coherent 
transport regime: 40 at any given energy, the quantum 
conductance is a constant value proportional to the num- 
ber of transmitting channels available for charge mobility, 
which are equal (in the absence of external leads) to the 
number of bands at the same energy. Therefore, the for- 
mation of effective dispersive bands in the investigated 
G4-wire is a signal of the intrinsic capability of the ma- 
terial of hosting electron " energy channels" available for 
charge migration, whithin continuous energy ranges. On 
the contrary, the energy spikes that would characterize 
a discrete spectrum would not be efficient for transport. 
The extended orbitals (see below) related to the effective 
bands turn out to be the corresponding viable "space 
channels" for carrier mobility, e.g. the pathways through 
which the carriers may migrate in the wire. 

A comparison between the manifolds of the isolated G- 
quartet and of the G4-wire highlights an increased den- 
sity of energy levels in the manifolds of the quadruple 



helix with respect to the tetrad. This effect is due to two 
factors: (i) the number of energy levels in each multi- 
plet; (ii) the in-plane and out-of-plane guanine-guanine 
and guanine-metal interactions, that couple the electron 
states stemming from the various structural elements (ei- 
ther from G or from K) in the supercell. The number 
of electron states in each manifold depends upon the 
number of molecules and ions in the unit cell: there- 
fore, whereas each manifold of the G-quartet contains 
four levels, the manifolds of the 3G4/K+ wire contain 
twelve levels from the twelve guanines, and additional 
twelve levels from the K + ions in the energy range where 
metal-base hybridization occurs. We come back to this 
point later when we discuss the DOS. The number of 
levels in a manifold is not the only element that deter- 
mines the details of the bandstructure: another key fea- 
ture is the band-width, which is instead controlled by 
the specific interactions. Whereas in the case of the G4 
tetrad only H-bonding plays a role, in the case of the 
3G4/K+ wire there are different contributions from H- 
bonding and stacking between the bases, and from the 
coupling between the bases and the metal ions. It is 
worth noting that the formation of dense manifolds and 
the consequent establishment of effective band-like po- 
tential conduction channels is a common characteristic 
of stacked H-bonded nucleobase aggregates. In fact, not 
only we identified the same features in guanine quadru- 
plexes both in the presence and absence of inner metal 
cations^ but similar energy manifolds were also detected 
in the electronic band-structure of two different DNA 
duplexes£& a i2£ For instance, in their simulation of an 
eleven base-pair poly(dG)-poly(dC)£ a sequence, de Pablo 
and coworkers found a HOMO-manifold derived from the 
eleven states of guanines. In that case, the topmost va- 
lence band had a bandwidth of only 40 meV, smaller than 
that calculated in our G4-wire (~ 700 meV). This con- 
firms that, while the splitting of the energy levels into 
multiplets is a fingerprint of stacked rotated H-bonded 
nucleobases, the peculiar characteristic of the manifolds 
depend on the intrinsic stacking properties (e.g. the 7r-7r 
coupling) of each type of helix. 

To investigate more closely the effects of the metal 
cations and of the metal-molecule interaction, we pro- 
jected^ the total density of states of 3G4/K+ (shaded 
gray area in Fig. 2b) on G (thin red curve) and K (thick 
blue curve) atomic orbitals. The G- and K-Projected 
DOS's (PDOS's) show that in a large portion of the spec- 
trum (below -2 eV and above 3.5 eV) only the guanine 
electron states contribute to the DOS. For instance, the 
LUMO peak at 3.5 eV has a pure guanine character, be- 
ing the convolution of the only twelve LUMO orbitals of 
twelve G's. The contribution of the K + ions to the total 
DOS is mainly localized at the top of the valence band 
(between -2 eV and eV). The presence of both guanine 
and potassium contributions to the DOS (as seen from 
the superposition of the blue and red peaks in Fig. 2b) 
at the top of the occupied bands is the consequence of 
a metal-molecule interaction. The HOMO peak of the 




FIG. 3: (a) Isosurface plot of the convolution of orbitals that 
represents the HOMO manifold of the G4-wire, side view, (b) 
Contour plot of the same orbital in a plane perpendicular to 
the tetrameric planes and containing the K + ions. 

3G4/K+ wire (see the double-peak indicated by an ar- 
row in Figure 2b) is due to a convolution of 24 electron 
states deriving from the coupling between the twelve G 
HOMO's and the twelve potassium orbitals. Indeed, one 
would expect that the K + ions would contribute to the 
DOS only with the filled p orbitals, after the s electrons 
are removed upon ionization of the system. However, if 
this were the case, by counting the levels in the manifolds 
we would find only nine completely occupied levels due to 
the metal (three p orbitals from each K ion), in addition 
to the guanine-based manifolds containing twelve levels 
each. Instead, the number of computed occupied levels 
(twelve more than those coming from the guanine count- 
ing, and all concentrated in the HOMO double-peak) and 
the occurrence of partial occupation for some of them 
(those around the computed Fermi level, see Figure 2b) 
indicates sp hybridization in potassium when inserted in 
the helical guanine complex. Therefore, potassium does 
not contribute separately with the 4s and the filled 3p 
shells, but with a partially occupied sp shell: 7 equiva- 
lent electrons in 4 orbitals, which become 6 electrons in 
4 orbitals when the system is ionized as the 3G4/K+ 
wire. Consequently, the highest-energy component of 
the HOMO double-peak has an occupation factor of 3/4 
the Fermi level is pinned at the top of the effective va- 
lence band (Fig. 2b) and the G4-wire in the presence of 
K + ions behaves as an intrinsically p-type doped system. 
The observation of delocalized holes at the topmost va- 
lence band is in agreement with the description of the 
electronic structure of potassium discussed above, where 
each ionized K + atom contributes to the whole electronic 
structure not with three fully occupied p-levels, but with 
four partially occupied sp-states. 

The spread sp-orbitals of potassium easily interact 
with the surrounding molecular orbitals of the tetrads, 
giving hybrid metal-molecule states. We visualized this 
effect by drawing the convolution of the effective HOMO 
of the wire. Figure 3a shows an isosurface of the HOMO 
charge density in the unit cell, obtained from the convo- 
lution of the 24 highest occupied orbitals that contribute 
to the HOMO double-peak. By comparison with Figure 
lb, it is possible to recognize a uniform distribution of 
the p orbitals deriving from the tetrads, but also the in- 
ner delocalized density stemming from the orbitals of the 
K + ions. By cutting the HOMO charge density with a 



plane perpendicular to the G-quartets, we further inspect 
the features of the metal- molecule interaction (Fig. 3b). 
Due to the high degree of charge derealization along 
the wire axis, the G4- wires may be described as good 
electron/hole channels for mobile charges. Two types 
of pathways for such mobile charges can be identified in 
Figure 3b and contribute to the conductivity channels. 
The first one stems from the low-energy component of 
the HOMO double-peak, is extended through the gua- 
nine core of the helix, and is due to the base-base inter- 
action. Similar channels are also observed for the other 
manifold-derived peaks (e.g. LUMO) and in the empty 
(K free) quadruple helix^i The second type of pathway 
is due to the high-energy component of the HOMO peak 
and results from the metal-base interaction. This hybrid 
pathway is centered around the potassium ions in the cen- 
tral cavity of the wire, and it clearly shows the coupling 
with the coordinated oxygen atoms of the G molecules. 
Hence, the inclusion of the metal cations inside the he- 
lix drastically influences the electronic properties of the 
system. The details of the metal-nucleobase interactions 
depend upon the nature of the cation included in the he- 
lix. In the present case of potassium, its coupling with 
the G bases modifies the topmost valence band, favoring 
the formation of extended orbitals along the stacking di- 
rection and enhancing the conduction properties of the 
empty guanine structured 

To get further insights into the metal-base coupling, we 
compared the 3G4/K+ wire with other two similar neu- 
tral systems: the empty quadruple helix (labeled 3G4) 
and the G4-wire with neutral K atoms (labeled 3G4/K). 
The structures 3G4 and 3G4/K do not describe real sys- 
tems, but are useful models for better understanding the 
3G4/K+ quadruplexes. The changes of the inner core 
do not basically modify the geometry (bond lengths and 
angles) of the quadruple helix, but affect the electronic 
structure. Figure 4 shows the comparison among the 
DOS (shaded gray areas) and the PDOS (thin red and 
thick blue curves) of the systems. To compare the three 
structures we aligned the bottom energy levels, and fixed 
the origin of the energy scale at the top of the valence 
band of the empty 3G4 helix (taken as reference), with 
the purpose to outline the effects of the inclusion of met- 
als into the empty guanine supramolecular structure. 

As already mentioned, in the DOS of all the studied 
G4-wires we observe the occurrence of broad G-derived 
peaks, reflecting the same guanine aggregation state. 
The empty tube has the valence band completely occu- 
pied and the Fermi level lays in the middle of the gap; 
in the presence of potassium (3G4/K and 3G4/K+), the 
Fermi level is pinned at the top of a partially occupied 
valence band. In the 3G4 case, the HOMO peak is the 
convolution of the HOMO's of guanines (Hq in Fig. 4), 
which are completely occupied; in the other two cases, in- 
stead, the HOMO derives from the interactions between 
the HOMO's of guanine (Hg) and the hybrid states of 
potassium (labeled K in Fig. 4), which are partially 
filled. One may consider the 3G4/K structure as ob- 
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FIG. 4: Total and Projected Density of States (DOS) of dif- 
ferent quadruple helical G4-wires: without inner metal (3G4) 
in the top panel, with K + ions (3G4/K + ) in the middle panel, 
with neutral K atoms (3G4/K) in the bottom panel. In each 
panel, the shaded gray area represents the total DOS, the 
thin red (thick blue) curve marks the contribution of gua- 
nine (potassium) to the DOS. Hg (Lg) labels the band-like 
peak deriving from the HOMO's (LUMO's) of single guanine 
molecules in the unit cell. The label K identifies peaks with a 
pronounced potassium contribution. The origin of the energy 
scale is consistently set at the HOMO of the empty column in 
the three plots. The 3G4 is an insulator with a large HOMO- 
LUMO gap. In the metal-containing quadruplexes 3G4/K + 
and 3G4/K, instead, the Fermi level is pinned at the top of 
the HG+K and of the K peaks, respectively, due to partially 
occupied electron states. 



tained by adding three electrons to the 3G4/K+. The 
inclusion of such electrons increases (from 3/4 to 7/8 fill- 
ing factor) but does not complete the occupation of the 
topmost band; thus, the helix with the atomic K would 
also be an intrinsic p-type doped wire. 

By analyzing the position and the width of the band- 
like peaks for the three systems, we can underline further 
details about the K-G interactions. As the K-derived 
states are found at the top of the valence band, far away 
from that range the spectra are identical, e.g. peaks 
a, b, e, f have the same width, shape and positions in 
the three plots. Also the G-like LUMO peak Lg remains 
basically unchanged. On the contrary, in the 3G4/K+ 
wire (central panel in Figure 4) the occurrence of a non- 
vanishing K-PDOS in the [-2;-l] eV range strongly mod- 
ifies the peaks c and d with respect to the empty wire: 
the electrostatic coupling between the s-like peak d and 



the sp orbitals of the K + ions shifts the position of this 
peak to lower energies and reduces its height. All the 
other peaks in Fig. 4 have a p-like character and are less 
influenced by the G-K + coupling. The displacement of 
peak d also modifies peak c, which is sharper and higher 
in 3G4/K + than in the free-standing tube. In the neu- 
tral 3G4/K helix, the electrostatic coupling between K 
and G is more effectively screened because of the larger 
number of electrons in the system: as a consequence, the 
shift of peak d is reduced with respect to the 3G4/K+ 
case, and peak d remains similar to that of the empty 
wire. Another difference that we remark is the compo- 
sition of the HOMO peaks: whereas for the 3G4 empty 
tube the HOMO is purely G-like (Hg), it is a mixed G- 
K-like double-peak (H G + K) for the charged 3G4/K+ 
wire. For the neutral 3G4/K wire, the HOMO (Fig. 4, 
middle panel) is split into two sharper peaks (Fig. 4, 
bottom panel): a lower energy component similar to the 
peak Hq and a higher component with a prevalent K 
character. The decrease in the degree of K-G orbital 
mixing m 3G4/K with respect to 3G4/K+ is most likely 
attributable to the increased occupation of the sp shell 
of potassium, which becomes more inert because more 
similar to a closed shell. 

These changes in the electronic structure due to dif- 
ferent charge states of the inner metals underline the 
delicate equilibrium that rules the stability and the 
mutual interactions in this hybrid organic/inorganic 
molecule/metal complex. Differently from M-DNA,13 
where the metal cation substitutes for an immino hydro- 
gen atom and is covalently bonded to a nucleobase, in G4- 
wires there is no direct charge sharing between the metal 
and the guanines. The stability and the electronic prop- 
erties of the system are ruled by the coordination ratio 
among the K + and the eight nearest-neighboring oxygen 
atoms, symmetrically located above and below the ion 
(Fig. lc). The metal cations stabilize the structure via 
electrostatic interactions, but also constitute an effective 
bridge between the electronic charge density around the 
oxygen atoms: the HOMO contour plot shown in Figure 
3b is the result of the 7r-7r coupling between the oxygens 
of consecutive tetrads, mediated by the electronic struc- 
ture of the K + ion. 

The possibility of changing the electronic properties 
of the G4-wire, by using different metal ions to stabilize 
the stack, is expected to be a powerful tool for tuning 
the conduction properties of the nanowires. We are cur- 
rently exploring the inclusion of other cations with dif- 
ferent outer-shell electronic configurations (e.g. the tran- 
sition metals), that may be exploited to further change 
the conductivity of G4-wires. 



IV. CONCLUSIONS 



The first principle study of the electronic and conduc- 
tion properties of infinite G4-wires, stabilized by K + in- 
ner cations, shows that the 7r-7r coupling among stacked 



planar tetrads is insufficient to induce the formation of 
dispersive bands in the wires. However, the presence of 
closely spaced energy levels leads to the formation of 
manifolds, whose density of states suggests a band-like 
behavior. The coupling among guaninc-localizcd molecu- 
lar orbitals, which may be easily induced by a weak exter- 
nal interaction, gives rise to extended electron channels, 
suitable to host mobile charge carriers along the wire. 

While the formation of split manifolds seems to be 
a general feature of base-base interactions in H-bonded 
stacked supramolecular nucleobase aggregates, the effects 
due to the presence of the metals depend on the identity 
of the cations. In the case of potassium, the inclusion 
of the cations enhances the conduction properties of the 
system, generating additional extended electronic chan- 
nels stemming from the metal-guanine interaction. The 
mixed sp-orbitals of K + hybridize with the HOMO's of 
guanine, and this coupling gives origin to a partially filled 
HOMO band, which makes the system equivalent to a p- 
doped wide-bandgap semiconductor. 

The above results, along with their attractive mechan- 



ical and self-assembly properties, suggest that G4-wires 
may be explored as viable DNA-based conductors for 
nanoscale molecular electronics. The computed proper- 
ties are the equivalent of those of a bulk material: of 
course, the actual behavior as a wire in a device setting 
depend on the specific device implementation, is affected 
from conditions such as electrode- wire and substrate- wire 
coupling, and leaves several open issues for further inves- 
tigation and exploitation. 
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